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Abstract. The Stackel transform is applied to the geodesic motion on Euclidean space, 
through the harmonic oscillator and Kepler-Coloumb potentials, in order to obtain ma- 
ximally superintegrable classical systems on TV-dimensional Riemannian spaces of noncon- 
stant curvature. By one hand, the harmonic oscillator potential leads to two families of 
superintegrable systems which are interpreted as an intrinsic Kepler-Coloumb system on 
a hyperbolic curved space and as the so-called Darboux III oscillator. On the other, the 
Kepler-Coloumb potential gives rise to an oscillator system on a spherical curved space as 
well as to the Taub-NUT oscillator. Their integrals of motion are explicitly given. The 
role of the (f lat/curved) Fradkin tensor and Laplace-Runge-Lcnz iV- vector for all of these 
Hamiltonians is highlighted throughout the paper. The corresponding quantum maximally 
superintegrable systems are also presented. 

Key words: coupling constant metamorphosis; integrable systems; curvature; harmonic oscil- 
lator; Kepler-Coulomb; Fradkin tensor; Laplace-Runge-Lenz vector; Taub-NUT; Darboux 
surfaces 
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1 Introduction 

The coupling constant metamorphosis or Stackel transform was formerly introduced in [H [2] and 
further developed and applied to several classical and quantum Hamiltonian systems in [3l HI El 
miT]. This approach has proven to be a useful tool in order to relate different (super)integrable 
systems together with their associated symmetries and to deduce new integrable Hamiltonian 
systems starting from known ones. 

For our purposes, the classical Stackel transform can be briefly summarized as follows OS]. 
Consider the conjugate coordinates and momenta q, p G with canonical Poisson bracket 
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{QijPj} = ^ij and the notation: 
N N 

1=1 1=1 

Let H be an "initial" Hamiltonian, Hjj an "intermediate" one and H the "final" system given 
by 

H = ^ + V{ci), Hu = ^ + U{cO, H=^ = ^ + V{ci), (1.1) 
/i(q) /i(q) [/ /i(q) 

such that fl = jiU and V = V/U . Then, each second-order integral of motion (symmetry) S 
of H leads to a new one S corresponding to H through an "intermediate" symmetry Sjj of Hjj. 
In particular, if S and Su are written as 

N 

S=Y. ^''i^)P^PJ + ^(q) = ^0 + VF(q), St/ = So + VFc/(q), (1.2) 

then one gets a second-order symmetry of H in the form 

S = So- WuH. (1.3) 

The aim of this paper is to apply the above procedure when the "initial" Hamiltonian H is 
the A^-dimensional (-/VD) free Euclidean motion, and when Hjj is either the isotropic harmonic 
oscillator or the Kepler-Coulomb (KC) Hamiltonian. It is well known that these three systems 
are maximally sup erintegr able (MS), that is, they are endowed with the maximum number of 
2N — 1 functionally independent integrals of motion (in these cases, all of them are quadratic 
in the momenta). These three systems and their MS property are briefly recalled in the next 
section, and we will see that the Stackel transform gives rise to several MS systems H that 
are defined on Riemannian spaces of nonconstant curvature. Moreover, we will show that the 
new potential V can be interpreted as either an (intrinsic) oscillator or a KC potential on the 
corresponding curved manifold. In this way, by starting from the Euclidean Fradkin tensor [8], 
formerly studied by Demkov in [9], and the Laplace~Runge-Lenz (LRL) A^- vector, the Stackel 
transform provides for each case its curved analogue (see [101 [TT| [12] and references therein) . 

In particular, we show in Section [3] that if Hjj is chosen to be the harmonic oscillator we 
obtain two different final MS Hamiltonians, for which H is endowed with a curved Fradkin 
tensor; these are a KC system on a hyperbolic space of nonconstant curvature and the so-called 
Darboux HI oscillator |13 ( I14 [ [15]. In Section |4] we take Hu as the (flat) KC Hamiltonian and 
the Stackel transform leads to other two different MS systems together with their curved LRL 
A^-vector; both of them are interpreted as intrinsic oscillators on curved Riemannian manifolds. 
Surprisingly enough, one of them is the A^D generalization of the Taub-NUT oscillator \16\ \T7\ 
[T8\ [T9\ [20l \2T\ [22} [23] . We stress that for some systems the dimension A^ = 2 is rather special as 
the underlying manifold remains flat, meanwhile for N >3 such systems are defined on proper 
curved spaces (see Sections 13.11 and 14. ip . This is similar to what happens in the classifications 
of 2D and 3D integrable systems on spaces of constant curvature (including the flat Euclidean 
one) [Ml [251 ESI [271 ISSl [291 EQl [31] which exhibit some differences according to the dimension 
and, in general, the 3D case is usually the cornerstone for the generalization of a given system 
to arbitrary dimension. 

As a byproduct of this construction, the "growth" of the Fradkin tensor and the LRL vector 
from their Euclidean "seeds" to their curved counterparts can be highlighted from a global per- 
spective. These results are comprised in Table [1] in the last section. Furthermore we also present 
in Table [2] the MS quantization for all of the above systems together with their "additional" 
quantum Fradkin/LRL symmetries. 



Superintegrable Oscillator and Kepler Systems on Spaces of Nonconstant Curvature 



3 



2 Harmonic oscillator and Kepler potentials on Euclidean space 

In order to fix a suitable common framework, we briefly recall the well-known basics of the 
superintegrability properties of the Hamiltonians describing free motion, harmonic oscillator 
and KC potentials on the A^D Euclidean space E^. 

As the "initial" Hamiltonian H (jl.ip for the Stackel procedure we consider the one defining 
the geodesic motion on plus a relevant constant a: 

H = ^p^ + a. (2.1) 

Obviously this is a MS system, and there are many possibilities to choose its integrals of motion. 
We shall make use, throughout the paper, of the following results. 

Proposition 1. 

(i) The Hamiltonian (12. ip is endowed with the following constants of motion. 

• (2N — 3) angular momentum integrals (m = 2, . . . , A^).- 

g(m) ^ ^ _ q.p.f^ S(^rn)= ^ {QiPj " QjPif , 

l<i<j<m N-m<i<j<N 

S(^) = S(;v) = L^, (2.2) 

where is the square of the total angular momentum. 

• iV^ integrals which are the "seeds" of the Fradkin tensor (i, j = 1, . . . , N): 

N 

Sij = PiPj such that Su = 2{H — a). (2-3) 

i=l 

• N integrals which are the "seeds" of the components of the LRL vector (i = 1, . . . , N): 

N N 

Si = ^Pk {qkPi - QiPk) such that ^ Sf = 2'L^{H - a). (2.4) 

k=l 1=1 

{ii) Each of the three sets {H, 3^""^, {H, 5(„,)} {m = 2,...,N) and {Su} {i = 1, . . . , N) is 
formed by N functionally independent functions in involution. 

(Hi) Both sets {H, S^"^\ ^(m), Su} and {H, S^'^\ S^^^-j, Si} (m = 2, . . . ,N and a fixed index i) 
are constituted by 2N — 1 functionally independent functions. 

As "intermediate" Hamiltonians Hjj (jl.ip we consider either the harmonic oscillator or the 
KC system. Since both of them are central potentials, the angular momentum integrals (12. 2p 
are valid for both cases, that is, = S'^'") and ^^^(m) = 5'(m) in (|1.2p . We recall that, 

in fact, the spherical symmetry of a central potential on E^ directly provides such (2A^ — 3) 
independent angular momentum integrals, so they characterize a quasi-MS system (321 [33] . 
However what makes rather special the harmonic oscillator and KC systems is the existence of 
one more independent integral, which is extracted from a new set of integrals that ensure their 
MS property and is related to the fact that these two systems are the only ones fulfilling the 
classical Bertand's theorem |34] . In this respect, each of the sets of integrals (12. 3p and ()2.4p 
gives rise to one known set of additional constants for the harmonic oscillator and KC system, 
respectively. 
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Proposition 2. 

(z) The harmonic oscillator Hamiltonian defined by 

Hu = + /3q=^ + 7 (2.5) 

has the {2N — 3) angular momentum integrals (|2.2|) together with N"^ additional ones given 
by the components of the ND Fradkin tensor {i,j = 1, . . . ,N): 

N 



Su,ij = PiPj + "^f^qiQj such that ^ Su,ii = 2{Hu - 7). 



i=l 



(ii) Each of the three sets {Hjj, S"*-™^}, {Hjj, 5(^)1 (?n = 2, . . . , N) and {Su,ii} (i = 1, . . . , A^) 
is formed by N functionally independent functions in involution. 

{Hi) The set {Hu, S^"^\ 'S'(m)) Su^a} {m = 2, . . . , N and a fixed index i) provides 2N — 1 func- 
tionally independent functions. 

Proposition 3. 

(i) The KC Hamiltonian given by 

^c/ = ^p' + A+e (2.6) 
2 |q| 

has the {2N — 3) angular momentum integrals (|2.2p together with the N components of the 
LRL vector {i = l,... ,N): 

N N 

Su,i = '^Pk{qkPi - QiPk) - such that S^j j = 21? {Hu -i) + 5^. 

k=l i=l 



{ii) Each of the two sets {Hu, S^"^^} and {Hu, S(^m)} {m = 2,...,N) is formed by N func- 
tionally independent functions in involution. 

{Hi) The set {Hu, S^^\ »S'(m), Sua} {m = 2, . . . , N and a fixed index i) is constituted by 2N — 1 
functionally independent functions. 

In the two next sections we apply the Stackel transform to each of these two MS systems. 
Notice that the proper isotropic harmonic oscillator arises whenever /3 = io'^/2 with frequency 00 
and 7 = 0, while the Kepler one corresponds to set 6 = —K and = 0. We remark that in 
this approach the constant a is essential in order to obtain a curved potential while the others 
(3, 7, 5 and ^ enter in both the kinetic and the potential term giving rise to MS oscillator/KC 
potentials on Riemannian spaces of nonconstant curvature, so that they can be regarded as 
classical "deformation parameters". 



3 Superintegrable systems from harmonic oscillator potential 

If we consider as the initial Hamiltonian H the free system (|2.ip and as the intermediate one 
the harmonic oscillator Hu ()2.5p . then we obtain the final Hamiltonian H 
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so that the relations (jl.ip read as 

^, = 2, V = a, C/ = 7 + /3q2, /i = 2(7 + /3q2) , V 



a 



7 + /3q2 ■ 

As far as the symmetries S = Sq + W (|1.2p are concerned, we find from Proposition [T] that 

^(m) = 0, So-ij = Sij, Wij = 0, 

while from Proposition [2] we obtain the elements Wjj for the decompositions of Sjj = Sq + Wjj, 
Wip^ = 0, Wu,i„,) = 0, Wu,^J = 2l3qiqj, 

where m = 2, . . . , N and i,j = 1, . . . ,N. 

Consequently, the Hamiltonian H (|3.ip is Stackel equivalent to the free Euclidean motion, 
through the harmonic oscillator potential, and its integrals of motion S come from (|1.3p and 
turn out to be 

S{m)^g(m)^ 5(^) = 5(„), S,,=p,pj-2f3q,q,H{q,p). (3.2) 

Thus we have obtained the (2A^ — 3) angular momentum integrals S'("^) and /S^^), together 
with N'^ ones, Sij, which form a curved Fradkin tensor. The quasi-MS property of H is ensured 
by the preservation of the (2A^ — 3) angular momentum integrals, that is, each of the two sets 
{H,S^'^^}, {H,S(^m)} = 2,...,A^) is formed by functionally independent functions in 
involution. Hence, from now on, we assume this fact and only pay attention to the additional 
constants Sij which characterize (j3.ip as a MS system. 

In order to perform a preliminary geometrical analysis of H we recall that, in general, any 
Hamiltonian of the form 

can be interpreted as describing a particle (with unit mass) on an AD spherically symmetric 
space Ai under the action of the central potential V(|q|) [H]. The metric and scalar curvature 
of M are given by |35] 



ds2 = /(|q|)2dq2, 

P 4)f (r)2 + /(r) {2f"{r) + 2(A - l)r-'f{r)) \ 

R = -^N-1}\ I, (3.3) 

where we have introduced the radial coordinate r = |q|. For general results on 2D and 3D 
(super)integrable systems on conformally flat spaces we refer to [36 1 137 1 [38]. 

Furthermore, the conformal factor /(|q|) = f{r) is directly related, under the following 
prescription, with the intrinsic KC and oscillator potentials on M: 

"Kc(r) := /' Uoir) := (3.4) 



that was introduced in [M] up to additive and multiplicative constants. 

With these ideas in mind, we now analyze the specific systems defined by H ()3.ip according to 
the values of the parameters j3 and 7. Notice that a is the constant which governs the potential, 
so to setting a = leads to geodesic motion on M, and that /5 must be always different from 
zero, since otherwise H is again the initial H. Therefore we are led to consider two different 
cases with generic a: (i) /3 / 0, 7 = 0; and (ii) /3 / 0, 7 / 0. 
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3.1 The case with /3 7^ and 7 = 0: a curved hyperbohc KC system 

If 7 = we scale H to deal with the Hamiltonian 

2 

Ukc = PH = ^ + 4- (3-5) 
Then, /(|q|) = |q| = so the metric and scalar curvature (|3.3p on M. reduces to 

= ii = - ^<^-;f (3.6) 

while the intrinsic KC and oscillator potentials (|3.4p on this curved space would be 
^Kc(r) = -^, Uo{r)=Ar\ 



The latter result shows that T-Lkg (|3.5p always defines an intrinsic KC potential on the space Ai. 
Nevertheless the curvature ()3.6p vanishes for N = 2, while the space is of nonconstant curvature 
for > 3. Therefore, for N = 2 the Hamiltonian must correspond to the usual KC system on the 
Euclidean space. This fact can be explicitly proven by applying to Hkc P-Sp the Kustaanheimo- 
Stiefel canonical transformation defined by |39l HQ] 

1/2 2\ ~ Pi9i-P292 - - P2qi+Piq2 „^ 

Qi = 1^{Qi- Q2), PI = 0,0 , Q2 = qiq2, P2 = 2 , „2 ' (3-7) 



pI+pI a l/~2 , ~2N , " 

^KC = ^,2^0, + -2—2 = obi +P2) + 



SO with canonical Poisson bracket {qi^Pj} = 5ij. In this way we recover the 2D KC Hamiltonian 

^I + pI 

2{ql+ql) ' ql+ql 2^^' ' ' 2^^^ 
and the five symmetries (|3.2p reduce to three integrals of motion, namely 
S'(^) = 5(2) =1? = {qiP2 - q2Pif = MqiP2 - hPif-. 

511 = Pi - 2qf'HKC = 2p2{q2Pi - qiP2) , „ - a, ^22 = -Sn - 2a, 

V^i + 92 

ctq2 

512 = S21 = P1P2 - 2gig2^KC = 2pi(giP2 - ^2^1) , „ „ ■ 

VQi + ^2 

Hence, by taking into account Proposition [3] for A^ = 2, we find that, under the above canonical 
transformation, the only angular momentum integral S^^^ is kept, while the four constants com- 
ing from the 2D Fradkin tensor reduce to the two components of the LRL vector: (Sn, 5*22) — "Si 
and (5i2, 521) §2- 

Consequently, a proper curved KC system arises whenever A^ > 3 and its full integrability 
properties can be summarized as follows. 

Proposition 4. 

(i) For N > 3, the Hamiltonian Hkc (|3.5p determines an intrinsic KC system on the hyper- 
bolic space of nonconstant curvature (j3.6p . 

{ii) Hkc is endowed with the {2N — 3) angular momentum integrals (12. 2p together with N"^ 
ones which are the components of an ND curved Fradkin tensor (i, j = 1, . . . , N): 

N 

Sij = PiPj — 2qiqj7iKC such that ^^5^ = — 2q and {Sii,Sjj}=0. 

1=1 

[Hi) The set {Hkc^ S(^m)j ^u} {m = 2,...,N and a fixed index i) is formed by 2N — 1 
functionally independent functions. 
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3.2 The case with f3 ^ and 7 7^ 0: the Darboux III oscillator 

If both /3, 7 7^ 0, we can write the Hamiltonian (j3.ip in the form 

Then the metric and scalar curvature (13.3P on the corresponding manifold M are given by 

9 / 9N 9 (Af-l)(2iV + 3(iV-2)Ar2) ^ ^ 

d^'=(l + W)dq^ R=-X^- ^^(^^^;,)3 ^ (3-9) 

and the intrinsic potentials (j3.4p read 



^Kc(r-) = -^^^^, ^o(r) = ^^. (3.10) 

In this way, we recover the A^D spherically symmetric generalization of the Darboux surface of 
type III [m HSl m] introduced in [HI [35]. Notice that the domain of r = |q| and the type 
of the underlying curved manifold depends on the sign of A [15] : 

A>0: i?(0) = -2AiV(iV- 1), rG[0,oo); 
A<0: i?(0) =2|A|iV(iV-l), re [0,1/^), 

where we have written the value of the scalar curvature (j3.9p at the origin r = 0. We stress 
that i?(0) coincides either with the scalar curvature of the A'^D hyperbolic space with negative 
constant sectional curvature equal to — 2A for A > 0, or with that corresponding to the A^D 
spherical space with sectional curvature equal to 2|A| for A < 0. 

By taking into account the above geometrical considerations and expressions (j3.10p . we find 
that T-L\ comprises both an intrinsic hyperbolic oscillator potential and a spherical one on M 
according to the sign of A. Strictly speaking the curved oscillator potentials arise by introducing 
the frequency oJ^ = — 2Aa and, in that form, the limit A — )• gives rise to the harmonic oscillator 
on E^, so A behaves as a classical deformation parameter governing the curvature and the 
potential. The MS property of Hx is then characterized by [T3| [T5]: 

Proposition 5. 

[i) The Hamiltonian T-Lx (|3.8p defines an intrinsic curved hyperbolic Darboux oscillator for 
A > and r € [0, 00) and a curved spherical Darboux one for A < and r G [0, 1/-^/|A[). 

{ii) Besides the {2N — 2>) angular momentum integrals (j2.2p . Hx Poisson-commutes with the N"^ 
components of the ND curved Fradkin tensor (i, j = 1, . . . , N) given by 

N 

Sij = PiPj — '^MiQj {T^x + ct) such that Sa = 2'Hx and {Su, Sjj} = 0. 

1=1 

(Hi) The set {Hx, S^^\ Sf^^^-j, Sa} (m = 2,...,N and a fixed index i) is formed by 2N — 1 
functionally independent functions. 



4 Superintegrable systems from the Kepler— Coulomb potential 

In this case we consider the initial Hamiltonian H ()2.1|) and the KG system Hjj (12. 6p for the 
intermediate one; this provides the final Hamiltonian H (jl.ip 

fr ^ lq|p^ , «lql (A 1 \ 
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where 



H = 2, V = a, U = — j — j — , /i = j — j , V — 



From Proposition [T] we obtain the decomposition of the symmetries S = Sq + W (I2.2p and ()2.4p 
{m = 2, . . . , N and i = l,...,N): 

while from Proposition [3] we find the one corresponding to Su = Sq + Wjj ()1.2p 

I ^1 

Therefore, the Hamiltonian H ()4.ip is Stackel equivalent to the free Euclidean motion, through 
the KC potential, and its integrals of motion S (ll.Sp are given by 

fe=l 

Hence, H ()4.ip is endowed with the (2A^ — 3) angular momentum integrals (j2.2p together 
with a curved LRL A^- vector with components Si. 

Notice that the parameter 6 cannot vanish in order to avoid the initial system H. Thus, 
similarly to the previous section, we study two systems covered by H (j4.ip : (i) 5 / 0, ^ = 0; 
and (ii) 5 / 0, ^ 7^ 0. 

4.1 The case with S ^ and ^ = 0: a curved spherical oscillator system 

If ^ = we have the Hamiltonian system 

= = ^|q|p^ + a|q|. (4.3) 

We stress that for A^ = 3 this system was early considered in [2]. The metric and scalar 
curvature (|3.3p give 

'dq^ ^,3(iV-l)(iV-2)_ 
|q| 4r 



and the intrinsic KC and oscillator potentials ()3.4p turn out to be 

UKcir) = Uo{r) = - 

y/r 4 

Hence T-Lq ()4.3p determines an intrinsic oscillator potential on M. However, for N = 2 the 
curvature is equal to zero, so this case should actually be the 2D harmonic oscillator. This can 
be proven by means of the canonical transformation 



?2 

11 = TTT?' Pi 



{piq2 - 2p2qi) + 92 - 91 + P2qi 



Q1+Q2-Q1) [VQi +Q2-Q1) 



l2={\/qf+qi-qi] , P2 



1/2 p2q2 - Pi [ Vli +Q2-Q1 



qf + qi- qi 



1/2 
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which is just the inverse of the canonical transformation (|3.7p : this yields the expected system 



The canonical transformation of the three symmetries (j4.2p gives 
5^^) = 5'(2) = = {q^p2 - q2Pif = ^{qiP2 - hPif , 

51 =P2{q2Pi - qm) + / T 2 " liPi-pf) + - q'i), 

VQi+12 * ^ 

52 = Pi{qiP2 - q2Pi) + , 'Ho = -P1P2 + aqiq2- 

V 9i + ^2 

Then the four components of the 2D Euclidean Fradkin tensor Sij are recovered, in the new 
canonical variables, from the set of constants (T-Lq, Si, S2) by setting 

Sn = 2{no + Si) =pI + 2aqi, S22 = 2(^o - ^1) =pI + 2aql 
S12 = S21 = 2S2 = P1P2 + 2aqiq2. 

Therefore the proper curved system arises whenever > 3, which yields the following 
Proposition 6. 

(i) For N > 3, the Hamiltonian T-Lq (|4.3p defines an intrinsic oscillator potential on the 
spherical space of nonconstant curvature (|4.4p . 

(ii) Ho Poisson-commutes with the {2N — 3) angular momentum integrals (j2.2p and with the 
components of the LRL N -vector given by (i = 1, . . . , N): 

N N 

Si = '^Pk {qkPi - qiPk) + pi^o such that ^ Sf = Uo- 2aV' . 
k=i i=l 

(Hi) The set {T-Lq, S^"^\ S(^m)^ Si} {m = 2,...,N and a fixed index i) is formed by 2N — 1 
functionally independent functions. 

4.2 The case with S and ^ ^ 0: the Taub-NUT osciUator 

We scale the Hamiltonian (j4.ip as 

nr, = CH= J'^^f'^ +^j^, rj = 6/(. (4.5) 
2(?7 + |q|) V+H 

Notice that the limit r/ — )■ reduces to the free Hamiltonian in Euclidean space. The metric and 
scalar curvature (j3.3p on the corresponding manifold M turn out to be 

d,^.!l±N,,., „(lM±£i^. (4.6) 

|q| /V ^ 4^(-^ + ' ^ ^ 

so that the domain of r = |q| in M depends on the sign of r]: 

r]>0: rG(0,oo); ?? < : r£[\r]\, 00). (4.7) 
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The intrinsic potentials (j3.4p are given by 



"KCM = -^/^, "oM = j^, (4.8) 

Consequently, T-L^^ defines two intrinsic oscillators, which are different systems according to (j4.7p . 

It is worth comparing g3]) with the Taub-NUT system [I6l[171[l8l[l9l[20l[2ll[22l[23] which 
can be written as [E^ 



rt Taub-NUT — , . 7]— jT + XT"; ^ I 1 + "j— T 

2(1 + 47n/|q|) 2(4m)^ \ |q| 



q|p2 ^ /i^|q|/(4m)^ ^ ^ ^V(4m) ^ 



2(4m + |q|) 2(4m + |q|) 2|q|(4m + |q|) 4m + |q 
The relationship with is established by setting 

^ = '"^' " = "2(4^' 
which gives 

_ |q|p^ _^ ^^|q|/(4m)^ ;uV(4m) 

rLn=dm + — TTT^ \ — I — TT + "^rr, \ — i — r 



(4m)2 2(4m + |q|) 2(4m + |q|) 4m + |q| 



so that we recover three terms in the "expanded" expression for ^Taub-NUT (|4.9p : namely, the 
kinetic term defining the geodesic motion on the Taub-NUT space (j4.6p . the insintric oscillator 
potential (j4.8p and the one which comes out by adding a constant to the oscillator potential. 
There is one missing term, the third one in (j4.9p . which corresponds to the Dirac monopole. 
However we notice that this can be derived from the angular momentum by introducing hyper- 
spherical coordinates in the form 



p'^ =Pr + r ^L^ and next + 

From this viewpoint, "H^ can be regarded as an A^D MS generalization of the Taub-NUT system 
which is recovered for 77 > 0, being the case with 77 < a different physical oscillator potential. 
The symmetry properties for Tirj are summarized in 

Proposition 7. 

(i) The Hamiltonian T-Lf^ ()4.5p characterizes two intrinsic oscillator potentials on the corre- 
sponding Riemannian space of nonconstant curvature (j4.6p according to (j4.7p . 

(ii) T-Lr} is endowed with the {2N — 3) angular momentum integrals l\2.2\i together with the 
components of the curved LRL N -vector given by {i = 1, . . . , N): 

N N 

Si = ^Pk{qkPi - QiPk) + 1] ri'^v such that Sf = 21? {Uj^ - a) + rj^Ul- 

k=l i=l 

{Hi) The set {Tir^, S^'^\ S^^^, Si} (m = 2,...,N and a fixed index i) is formed by 2N — 1 
functionally independent functions. 
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Table 1. Maximally superintegrable classical oscillator and KC Hamiltonians in N dimensions. 



• Geodesic motion on Euclidean space 
* Common (2A'^ — 3) angular momentum integrals of motion 

l<i<j<m N -m<i<j<N 

* "Seeds" of the A'^D Fradkin tensor * "Seeds" of the LRL A''-vector 



Sij — PiPj 



Si ^^^Pk {qkPi 



q%Pk} 



• Harmonic oscillator 



* Flat A''D Fradkin tensor 
Su,ij = PiPi + "^Pqiqj 



• Euclidean KC 



2^ +R + ^ 



* Flat LRL iV-vector 

AT 

Su,i ^'^''P^ ^ ^ TT 

fe=i 



• Curved hyperbolic KC {N > 3) 

2q"^ q^ 

* Curved A'^D Fradkin tensor 

Sij = PiP] - 2qiqjHKC 



• Curved spherical oscillator {N > 3) 
Wo = ^ |q|p^ +a|q| 

* Curved LRL A'^- vector 

AT 

Si = '^Pk {qkPi - qiPk) + y-Mo 



Darboux III oscillator 

p^ \oi q^ 



2(1 + Aq2) 1 + Aq2 
* Curved A''D Fradkin tensor 



Sij = PiPj - 2\q,qj (Hx + a) 



• Taub-NUT oscillator 

,2 



|q|p 



+ 



Q|q| 



2(?? + |q|) ^+|q 



* Curved LRL A''- vector 

JV 

Si = '^Pk [qkPi ~ q%Pk) + J? Pi 'H,, 
fc=i l*^' 



5 Outlook and superintegrable quantization 

So far we have obtained and interpreted four MS classical Hamiltonian systems on Riemannian 
spaces of nonconstant curvature by starting from free motion on E^^ and applying the Stackel 
transform through the harmonic oscillator and KC potentials. The main results here obtained 
are displayed in Table [1] where the transition from the "seeds" of the Fradkin tensor and the LRL 
vector up to their curved analogues is laid bare by reading the table through its two columns. 
Recall, however, that the Darboux III and the Taub-NUT oscillators give rise, each of them, to 
two different physical systems according to the sign of the parameters A and r/, respectively. 

Some related comments are in order. All the Hamiltonians shown in Table [T] are constructed 
on spherically symmetric spaces so that they are endowed with an 5o{N) Lie-Poisson symmetry. 
In particular, let us consider the generators of rotations Jij = q-ipj — qjpi with i < j and 
i,j = 1, . . . , iV which span the 5o{N) Lie-Poisson algebra 

{Jij, Jik} — Jjkj {Jij J Jjk} — Jikj {Jikj Jjk} — '^ijj i ^ j ^ k. 

Then the "common" (2N — 3) angular momentum integrals S*-™^ and (j2.2p can be written 
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Table 2. Maximally superintegrable quantum oscillator and KG Hamiltonians in N dimensions. 



* Common {2N 

l<z<j <m 



3) quantum angular momentum operators 

N — m<i<j<N 



• Quantum hyperbolic KC (A'^ > 3) 


• Quantum spherical oscillator (A'' > 3) 


1 ~ 2 <^ 


= i |q|p^ + Q|q| 


* Quantum 7VD Fradkin tensor 


* Quantum LRL Af- vector 


Sij = PiPj ~ 2q^qj'HKC 


1 ^ 1 ^ - 
5« = - ^Pfe {qkPi - qiPk) + 2 5Z ^^'^P^ ~ iiPk)Pk + ■Pm'^o 
fc=i fe=i l*^' 


N 

Y,S^^ = -2a 


^ 1 


1=1 


1=1 



• Quantum Darboux III oscillator 

il ^ 1 -2 Aaq^ 
^ " 2(1 + Aq2) ^ 1 + Aq2 

* Quantum A'^D Fradkin tensor 



• Quantum Taub-NUT oscillator 



|q| .p2 4^ , "I'll 



2(^7 +|ql) r7+|q| 
* Quantum LRL A'^- vector 



Sij = PiPj - 2\qiqj ['H\ + a 



Sr, = 2?i> 



1 1 

Si = - ^Pk {qkPi - QiPk) + 2 X] ^^''P' ^ qiPk)Pk + rj rrvTiv 

fc=i fe=i 1*1 

^ 1 

^ S? = 2L2(H^ - a) + 77^?^^ + -{N ~ ifh^ii^ - a) 



as the quadratic Casimirs of some rotation subalgebras so(m) C so(A^): 

l<i<i<"i N-m<i<j<N 

with 5^^-* = >S'(7v) = being the quadratic Casimir of so(A'^). In this respect, we also notice 
that all the LRL constants of motion Su^i, Si) given in Table [T] are transformed as A^- vectors 
under the action of the generators of so(A^) (as it should be): 

{Jij, Sk} = SikSj — 6jkSi. 

Furthermore, all of these systems possess an s[(2,M) coalgebra symmetry as well [U]. If we 
denote J_ = , J+ = and J3 = q • p we have that 

{J3,J+} = 2J+, {J3,J_} = -2J_, {J_,J+} = 4J3, 

and the common integrals ()2.2p are just the mth (left and right) coproducts of the Casimir of 
s[(2,]R). This set of (2A^ — 3) integrals is "universal" for any Hamiltonian function defined by 
H = ^(q^,p^,q • p) so that this always provides, at least, a quasi-MS system \32\ I33j . There- 
fore the Hamiltonians shown in Table [1] are distinguished systems since they have "additional" 
symmetries. 

To end with, we shall present the MS quantization of the four curved classical systems. For 
this purpose, we remark that the MS quantization of the Darboux III oscillator has been recently 
obtained in [45], and the corresponding quantum dynamics has been fully solved for A > (the 
case with A < is still an open problem). This quantization has been obtained by applying 
the so called "Schrodinger quantization" procedure [l6], and its relationship with the Laplace- 
Beltrami and position-dependent-mass quantizations has been established in j47] by means of 
similarity transformations (see also |48j). 
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Therefore, let us consider the quantum position and momenta operators, q, p, with canonical 
Lie bracket = ihSij. The resulting MS quantum Hamiltonians are summarized in the 

following final statement. 

Proposition 8. Let Ti he one of the classical Hamiltonians given in Propositions^^ 

[i) The Schrodinger quantization ofH and its quantum, symmetries are given in Tabled 

(a) The quantum Hamiltonian % is endowed with {2N — 3) quantum angular momentum ope- 
rators and Si^rn), such that {Ti,^^™)} or {U^Si^rn)] {m = 2, . . . ,N) is a set of N 
algebraically independent commuting ohservables. 

{Hi) If % a (HkCi^a) then it commutes with the N"^ components, Sij, of a quantum Fradkin 
tensor (i, j = 1, . . . , A^). The set {7^, S^"^\ 5'(m), Su} {m = 2, . . . , N and a fixed index i) is 
formed by 2N — 1 algebraically independent commuting observables. 

(iv) When % € {'Ho,'Hri) this commutes with the N components, Si, of a quantum LRL vector 
(i = 1, . . . , N). The set {T-L, 5^™-*, S(^„^^, Si} {m = 2, . . . , N and a fixed indexi) is constituted 
by 2N — 1 algebraically independent commuting ohservables. 
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